Spatial adiabatic passage in a realistic triple well structure 
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We investigate the evolution of an electron undergoing Coherent Tunnelling via Adiabatic Passage 
(CTAP) using the solution of the one-dimensional Schrodinger equation in both space and time 
for a triple-well potential. We find the eigenspectrum and complete time evolution for a range 
of different pulsing schemes. This also provides an example of a system which can be described 
with the tools from both quantum optics and condensed matter. We find that while the quantum 
optics description of the process captures most of the key physics, there are important effects which 
can only be correctly described by a more complete representation. This is an important point 
for applications such as quantum information processing or quantum control where it is common 
practice to use a reduced state space formulation of the quantum system in question. 

PACS numbers: 



I. INTRODUCTION 



The fabrication and control of mcsoscopic quantum 
systems is of continued interest from both a fundamental 
and practical point of view. One of the most exciting as- 
pects of this is the convergence between quantum optics 
and quantum electronics. While these have been tradi- 
tionally treated as separate but related fields, the raft 
of new experimental and theoretical results have shown 
many important links between these topics. This is espe- 
cially true in quantum information processing and quan- 
tum control where the techniques and concepts used in 
one physical system are regularly applied to the other. 
An important aspect of this convergence of quantum op- 
tics and electronics is the interoperability of the tools 
used to describe their respective systems. 

To explore the convergence between quantum optics 
and quantum electronics, we consider the recent analy- 
sis of STIRAP (Stimulated Raman Adiabatic Passage), 
which was originally developed in the quantum optics 
framework^, applied to spatial coherence of electrons, in 
which case it is termed CTAP (Coherent Tunnelling via 
Adiabatic Passage)^. This is an example of a non-trivial 
time dependent problem which has both a quantum op- 
tics and condensed matter version and which is of inter- 
est in both communities. The CTAP protocol has been 
suggested as a means to transport electrons from one 
potential well to another via adiabatic manipulation of 
the ground state wavefunction. The transport of spins 
using this protocol has been proposed as the basis of 



quantum computing schemes^i^. Similarly, transporting 
single atoms^ and Bose-Einstein condensatesii^ within 
harmonic traps have also been investigated. CTAP has 
recently been demonstrated using photons in triple-core 
optical waveguides, using modulated separations between 
the cores2ii£iiii. 

The development of the CTAP protocol was carried 
out directly from the quantum optics framework which 
translates to an ideal localization assumption. In this 
work we relax this assumption and investigate the pro- 
tocol in a manner more appropriate to a realistic quan- 
tum nanoelectronic/solid-state setting, and highlight the 
differences and commonalities of the approaches to each 
physical context. The results of our analysis are directly 
applicable to the implementation in the quantum dot set- 
ting, including gate define d^^i^^'^'^i^^'^^i^'^'^^ and donor 
based Coloumb confined^ single electron systems. 

To this end, this paper necessarily focuses on describ- 
ing few state quantum systems using tools common to 
either the quantum optics or condensed matter commu- 
nities. The theoretical treatment of CTAP and other 
related scheme o^^i^°'^-^i^^'^'^i^^i^^ that involve the coher- 
ent behaviour of electrons often employs the finite state 
approximation (FS), also variously known as the tight- 
binding or modal approximation. As such, the electron is 
located at one of a number of localised sites and can tun- 
nel coherently between these sites. The evolution of the 
system is then solved using matrix mechanics using the 
finite state space. In this model, the mathematics of the 
system evolution is often equivalent to a few state quan- 
tum optics problem, especially when considering donor 
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electrons in semiconductors or quantum dots containing 
few electrons. This approach has become particularly 
important given recent work on quantum control and 
quantum computing in both quantum optics and solid- 
state^i^. 

We investigate the generic problem of the link between 
the finite state model parameters with those that can 
be determined from conventional solid-state physics and 
the solution to the Schrodinger equation in one dimen- 
sion using wave mechanics. We use CTAP as a specific 
example where the appropriate time variation of the po- 
tentials is calculated for a series of square wells separated 
by finite barriers. We discuss the analysis of the result- 
ing eigenspectrum and then perform ID time dependent 
calculations to investigate the system evolution. 

At this point, it should be noted that the solution of 
the spatially varying wavefunction for an arbitrary 3D 
potential is the more sophisticated and ultimately cor- 
rect treatment. At best, the FS treatment is a physi- 
cally motivated approximation for treating the shuttling 
of electrons between quantum dots or donor sites. Here 
we choose a mid point between FS and full 3D, where 
the dynamics are dominated by the potential variation in 
only one of the three dimensions, the other two dimen- 
sions consisting of tight confining potentials. Most of the 
essential physics involving barrier penetration and finite 
confinement enter the problem at this ID level. Further 
effects introduced in the full 3D problem are mostly sys- 
tem specific, depending on both the characteristics of the 
host lattice and the nature of the confining potentials. 
Investigating the similarities and differences between the 
Schrodinger wave (SW) description, even in ID, and the 
finite state approximation (FS) provides important in- 
sight into the validity of this approximation. 

While this paper focuses on solutions to the one dimen- 
sional Schrodinger equation, the analysis and methodol- 
ogy applies to many real systems, as coherent tunnelling 
can often be treated as a one dimensional problem. In 
GaAs quantum dots for instance, the separability of the 
spatial degrees of the wavefunction often allows the dy- 
namics in one dimension to be treated independently of 
the confinement of the wavefunction in the other two di- 
mensionsH. A notable exception is the charge states as- 
sociated with isolated dopant atoms in semiconductors, 
where the Coulombic nature of the trapping potential in- 
troduces divergences which must be handled with some 
care. In this case, our analysis still provides some insight 
as similar qualitative behaviour is seen in these systems 
when compared to other forms of trapping potentials. 

In the following section we discuss the links between 
the approaches used in quantum optics (the FS model) 
and those commonly used in condensed matter physics 
(the SW model). We then introduce CTAP in section Hill 
as an example of a nontrivial adiabatic process which can 
be modelled using either formalism. As CTAP is an adi- 
abatic process, much insight can be gained by plotting 
the eigenspectrum. In section IIVI we find the eigcnen- 
ergies and wavcfunctions using the SW formalism and 



use this to estimate the adiabatic time scale. Using the 
FS description of CTAP, the population of the central 
well is predicted to be exactly zero in the adiabatic limit. 
Using the SW formalism, section |Vl we can put quanti- 
tative bounds on this population. Finally in section [Vll 
full time dependent evolution is calculated for the CTAP 
process and this is used to investigate the fidelity of trans- 
fer as a function of switching speed. This allows direct 
comparison between the FS and SW descriptions. 

II. THE DOUBLE WELL SYSTEM AND THE 
LINK TO QUANTUM OPTICS 

The canonical system for investigating two-state quan- 
tum mechanics in the solid-state is the double well po- 
tential, whereas in atom optics it is the two-level atom 
in the rotating wave approximation. Much of the mathe- 
matical properties of these two descriptions arc trivially 
equivalent, which leads to them being often used almost 
interchangeably. In this paper, we are primarily con- 
cerned with the points where this equivalence becomes 
non-trivial and ultimately breaks down, especially in the 
case of multi-well systems. To motivate the discussion 
and introduce the relevant formalism, we will first point 
out several important features of the two-state and/or 
double well system as background. 

Consider a one dimensional system consisting of two 
potential wells separated by a finite potential bar- 
rier, Fig[T] We assume that the potential barrier between 
the wells is sufficiently high that each well contains at 
least one 'bound' state. We also assume that the wells 
are narrow enough that the higher lying eigenstates of 
the system are sufficiently removed from the lowest two 
states that we can approximate the system as a two-state 
system. The Hamiltonian governing the evolution of the 
lowest two eigenstates can be written as 

in a basis (|L) and \R)) comprised of the ground states 
of the left- and right-hand wells respectively. The cou- 
pling between states k is just a function of the height and 
width of the finite barrier and El and correspond to 
the energies of each well defined separately, given an in- 
finite barrier. It is at this point that the analogy with a 
two-level atom in the rotating frame from quantum op- 
tics is made. In a driven two level atom problem, the 
off-diagonal terms arc given by the intensity of the driv- 
ing field (the Rabi frequency) while the diagonal terms 
come from the energy mismatch between the bare tran- 
sition and the frequency of the driving field (the detun- 
ing). In the square well problem, the link between this 
reduced state space and the original spatially dependent 
wavefunction is then contained in the calculation of k for 
a finite barrier. 

The solution of the double well problem is most of- 
ten constructed by matching boundary conditions of 
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which allows us to construct the Hamiltonian of the sys- 
tem on a discrete lattice, 



h w h 

FIG. 1: The canonical two square well, finite barrier problem. 




FIG. 2: The energy gap between the ground and excited states 
(AsAs) as a function of the barrier height Vo in units of E* , for 
a double well system comprised of wells of width L separated 
by a barrier of width w, as shown in Fig. [T] 



the solution to the Schrodinger equation in each region 
separately^. In the case of the double well, the solution 
does not have a closed form but can be evaluated numer- 
ically. The calculation of k then reduces to finding the 
energy gap (AgAs) between the two lowest (symmetric 
and antisymmetric) eigenstates of the double well sys- 
tem. Fig. [5] gives the value of Asas for two wells of 
width L separated by a barrier of width w and height Vq. 
The energy is expressed in normalized units {E*) which 
corresponds to the ground state energy of an infinite well 
of width L. 

A second method of solving the ID Schrodinger in- 
volves direct numerical discretization of the wavefunc- 
tion, and will be used later for solving an arbitrary ID 
potential problem. If we discretize the region of interest 
into n points xi,X2, ■ ■ ■ separated by Aa;, the wave- 
function can then be computed at each point -0^ = 'ip{xk)- 
We then express the second derivative as a finite differ- 
ence equation 



Hswixk) = V{xk) 



tpk+i - 2V'fc -I- ipk-i 
(Ax)2 



(3) 



and express the Schrodinger equation as a matrix equa- 
tion. For our purposes we truncate the matrix at ei- 
ther end, which corresponds to infinite potential bound- 
ary conditions, V{xo) = V^Xn+i) = oo. Diagonalizing 
this Hamiltonian matrix gives the eigenvalues, the low- 
est ones of which are good approximations to the system 
eigenenergies. Similarly, the wavefunctions will be the 
eigenvectors corresponding to these eigenvalues. While 
this method is not particularly sophisticated, it does al- 
low the inclusion of arbitrary potential shapes which will 
be important later. In addition, because the resulting 
matrix is tridiagonal, this diagonalization process can 
be performed quickly and efficiently with modern sparse 
solving routines^^. 



III. CTAP IN A REALISTIC POTENTIAL 

To illustrate the link between the finite state and 
Schrodinger wave approaches, we take as an example a re- 
cently proposed technique for moving electrons spatially 
using coherent, adiabatic evolution. As this scheme is the 
direct electrical analogue of STIRAP, the existing anal- 
yses specifically uses the FS model. Here we investigate 
the state evolution and pulse design for this system using 
the SW approximation. The full treatment of the spatial 
extent of the electron is ultimately the more correct de- 
scription as it does not assume that the electron is tightly 
bound to one of the sites included in the finite state basis 
set. 

We begin by reviewing the FS treatment of CTAP. 
Given a single electron which can be localised in one of 
three distinct spatial regions, we describe this in a basis of 
three states which are nominally the ground state of the 
electron at each of the three sites. These spatial locations 
are analogous to energy levels in a multi-level atom in 
quantum optics. Obviously the first approximation here 
is that the excited states of the quantum wells can be 
ignored in this process. This means any physics derived 
from the population of these higher levels is essentially 
lost, but we will discuss this later. 

We assume that the tunnel barrier between the dots 
can be modulated in such a way that we may vary the 
coherent tunneling rates (or matrix elements) between 
dots. This process is directly comparable to laser medi- 
ated transitions between levels (in the rotating frame) of 
a multi- level atom (see Fig. ^ . In the finite state model, 
the Hamiltonian is given by 
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(Ax) 



(2) 
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FIG. 3: Correspondence of the (a) quantum optics description 
of a three level atom with laser transitions which is equivalent 
to the finite state (FS) approximation and (b) a single electron 
undergoing coherent tunnelling between one of three finite 
wells (SW). The tunnel matrix elements between states are 
Ki2 and K23, while the single-photon detuning (central well 
offset energy) is given by A. 



where A is the offset energy (detuning) of the central site 
relative to the end sites, K12 and K23 are tunnel matrix 
elements and wc have set the energies of the first and last 
sites equal. 

The eigenstates of this Hamiltonian arei*^ 

|2?+) = sinOi sin 02|a) + cos 82I6) + cos 81 sin 02|c), 
= sin6i cos 92|a) — sin92|6) + cos 61 cos 62IC), 
\Vn) = cosGila) +0|5) -sinGilc), (5) 

where 



61 = arctan (K12/K23) , 
1 



02 = — arctan 



V(2ki2)2 + (2«;23)2 /A . (6) 



The eigenenergies of these states are 
A 1 



£± = -±-v/(2ki2)2 



(2^23)2 +A2, 



£0 



0. 



(7) 



Applying a specific pulse sequence to modulate the co- 
herent tunnelling allows the transfer of the electron be- 
tween the two extrcma wells. If the barrier closest to the 
electron is lowered first, the system adiabatically follows 
the ground state eigenstate spending some time in the 
central well. This we refer to as the intuitive direction 
(ID). Reversing the pulse sequence is termed the counter- 
inuitive direction (CID) and similarly results in a transfer 
of the electron but this time via the 1st excited state of 
the system. This can be seen directly from the eigen- 
states Eq. where the CID path transforms from \a) 
at t = to |c) at t = iinax and hence the goal of the CID 
sequence is to maintain the system in |I?o) throughout 
the evolution. 

The original work on CTAP— assumed that the tun- 
nelling rates between wells could be arbitrarily con- 
trolled, whereas here we are specifically concerned with 
how to modulate the barrier heights to achieve an ar- 
bitrary tunnelling rate variation. It is by investigating 



this link between the control of the barrier height and 
the tunnelling rates and associated issues that we can in- 
vestigate the link between the FS and SW descriptions. 
In this case, the SW version of the problem is described 
by a ID potential consisting of three wells, separated by 
finite barriers (i.e. Fig. [HJd). The discrete form of the 
Hamiltonian of the system is then given by Eq. ([5]). 

If we start with the conventional CTAP pulsing 
scheme, the coupling between states is varied in time in 
a Gaussian fashion, 

U\ / [t-(tmax+^)/2]2\ 
Ki2(t) = Kniaxexp<^ ^ ^ (8) 

K23{t) = Kmaxexp <^ ^ '> , (9) 

where K,nax is the maximum coupling, imax is the pulse 
time and a controls the width of the pulses in time. 
To achieve an equivalent modulation using finite wells, 
we need to vary the barrier heights starting from some 
maximum barrier height (Vmax), down to some minimum 
height (Vjnin) and back up to the maximum. The CTAP 
pulses are then constructed by applying identical pulses 
to both barriers with some time delay, cr, between them. 

To find an approximate barrier modulation pulse to 
replicate the behaviour of the Gaussian pulses, Eq. ([S]), 
we can consider the conventional result that the proba- 
bility of tunnelling through a finite (but large) barrier of 
height Vq and width w is proportional to expl—wy/Vo]. If 
we naively equate this to the coupling energy we obtain 



exp[ 



Vol (X exp <^ - 



[t - (tn,ax±g)/2]^ 



(10) 



and normalizing for the minimum and maximum barrier 
heights gives 



Vit) 



(^max "l~ Cr)^ 



(2t-<,„axT^T) +Knin, (H) 



which we designate as pulsing scheme I. Fig. UJa) shows 
the barrier heights given by Eq. pT|) as a fraction of the 
pulse time for a = tmax/16, T4iin = 10 and Knax = 1000. 

To derive a more rigorous pulsing scheme (ll), we can 
take the previous calculation of coupling as a function 
of barrier height (Fig. [2]) and use interpolation to calcu- 
late the required barrier height for a given coupling k. 
Fig. m^b) shows this pulsing scheme for a barrier width 
oi w = 0.2L. The functional form of this pulse is very 
similar to that of pulse I but with a smooth roll-off at the 
pulse maximum. While pulse II is more difficult to cal- 
culate (as it requires the numerical solution to the finite 
well problem) this calculation need only be performed 
once with the pulses being then constructed via lookup 
tables and interpolation. 

The pulses considered so far start and finish with the 
barriers both 'high'. It is also possible to construct a 
pulse where only one barrier is initially high, the other 
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FIG. 4: Four diflerent pulsing scliemes for modulation of a 
barrier height to result in a CTAP transfer between wells. In 
each case, the solid (dashed) curve is the height of barrier 1 
(2) as a function of pulse time. Note the log scale for pulse IV. 



barrier (which is not required to confine the electron ini- 
tially) can be started in its 'low' stato^. An example of 
this is 

Viit) = (K.ax-Kiin)cOs2 ( ) +K.i„ (12) 

V2it) = (Kiax-KihOsin^ ( J-^) (13) 

\ ^ ''max / 

which we designate pulse III, see Fig.|4l[c). In this pulse 
scheme there no longer exists an ID pulse and the CID 
pulse corresponds to following the 2nd excited state of 
the system. In a system in which the native barrier po- 
sition is low and voltages need to be applied to raise the 
barrier height, an alternating pulse of this type may be 
simpler to implement. It also does not require timing de- 
lays between the pulses, as the pulse overlap stems from 
the time reversal symmetry of the pulses. 

While pulse III is relatively simple to describe and/or 
implement, it is not optimal in the sense that the great- 
est variation in barrier height is exactly where there is 
greatest overlap, resulting in a very narrow energy gap 
for adiabatic evolution. A better pulsing scheme would 
be to construct pulses which display a trigonometric vari- 
ation in coupling rather than barrier height. Again using 
Fig. [2] to obtain such a pulse (IV) we obtain the pulse 
scheme shown in Fig. IH^d), which requires much finer 
control in the initial and final stages of the pulse but re- 
sults in a much smoother variation in the coupling over 
time and large adiabatic gap. 



IV. EIGENSPECTRA AND ADIABATIC 
CRITERIA 

Using the method of section [III we can obtain the 
eigenenergies and eigenstates of the triple well system 
as a function of pulse time. In this way we can construct 
the eigenspectrum, which provides information about the 
adiabaticity of the system, as well as determining the 
population distribution at each point in time. While in 
principle the FS formalism provides this same informa- 
tion, using the SW method allows the inclusion of real- 
istic spatially varying potentials. This, in turn, allows 
calculation of the higher lying states and the spatial dis- 
tribution of the wavefunction itself. 

In Fig. [5] we show the eigenspectrum of the finite well 
system (consisting of 3 wells of width L separated by 
two barriers of width w = 0.2L) undergoing evolution 
due to modulation of the barrier heights according to 
pulse II. In this plot, the manifolds of three states are 
visible at progressively higher energies. Each of these 
manifolds corresponds to a bound state in a well of width 
L. The width of the wells has been chosen so that there 
is a large separation between each group of three states 
which guarantees that the adiabatic evolution will leave 
the particle in the ground state of the occupied well at 
the end of the pulse sequence. 

For this and subsequent calculations involving pulses I 
& II, the bottom of the central potential has been raised 
such that V{-L/2 <x< L/2) = Q.IE*, to separate the 
energies of the states within each manifold. The sign of 
this offset does not qualitatively change the following re- 
sults, provided that the energy of the lowest 3 eigenstates 
does not approach those of the next highest manifold of 
states. 

For pulses III & IV, the central well offset is unneces- 
sary as the initial condition (with only one barrier high) 
naturally splits the degeneracy of the first three eigen- 
states. In contrast, for all pulses sequences the left and 
right wells are identical, resulting in no offset between 
their ground state energies. This has been shown to be 
an important rcquirementii^i^ as variation in energy be- 
tween the initial and final well results in crossings in the 
eigenspectrum, which limits the adiabatic evolution. 

Fig. [5] shows the eigenstates of the three well system 
at three different points during pulse sequence II, whose 
eigenspectrum is shown in Fig. O The eigenstate corre- 
sponding to the CID pathway for CTAP is given by the 
1st excited state (the solid line) in each subfigure. Note 
that the electron is well localised at the beginning and 
end of the pulse, but the spatial location of the eigen- 
states have swapped. At the halfway point of the pulse 
sequence, the eigenstates are superpositions of all three 
well states as well as having some non-negligible popu- 
lation within the barriers. In contrast to the FS model 
the population in well 2 is not zero, even in the limits of 
long time and ideal coupling, due to the extended nature 
of the wavefunction. However in common with the FS 
model, the central well population is strongly suppressed 
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Time [t/tmax] 

FIG. 5: Eigenspectrum of the three well system as a function 
of time within pulse sequence II. The insert shows the three 
lowest eigenvalues in the central region of interest where the 
characteristics of the eigenstates are changing fastest. 
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FIG. 7: Eigenspectrum of the three well system as a function 
of time for (a) pulse III and (b) pulse IV. Note the very small 
energy gap for pulse III, due to the smaller tunnelling rate 
at the midpoint of the pulse sequence. 
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FIG. 6: Eigenstates of the three well system as a function of 
time within pulse I plotted at (a) t/tmax = 0.2, (b) t/tmax = 
0.5 and (c) t/tiaax = 0.8. The eigenstates are the ground state 
(dashed), 1st excited (solid) and 2nd excited state (dotted) 
respectively, where the CTAP pulse corresponds to the 1st 
excited state. Note, the potential is shown in arbitrary units 
for comparison. 



for the CID pathway state. We discuss this in more detail 
in section IVl 

We can perform the same analysis to obtain the eigen- 
spectrum of the other pulse schemes, which are shown in 
Fig. [71 Note that, given the similarity between pulse I 
and II; the eigenspectrum from Fig. [5] is not presented 



here. From Fig. [7] we immediately see that pulse III 
is indeed a poor choice as it results in a very small (but 
non-zero) gap between the first and second excited states. 
In contrast, pulse IV provides a smooth variation in the 
energies of the system with a similar gap size to that 
resulting from pulses I and II. 

Once we can calculate the eigenspectrum of the system, 
we can also estimate a time scale over which the system 
evolution is adiabatic. To estimate this, we use 



{m\'-i\m) 



\{d^,{t)\H\ct^,{t)) - mmmw 

(14) 

as the adiabaticity of a transfer between eigenstates 4>i 
and 0/^. We can then define the adiabaticity of the 
CTAP process as 

^CTAp(i) = max{^[0i(t),02(t)],^[<^2(t),03(t)]} (15) 

where the (/)'s are the first three eigenstates of the system 
(ordered by energy) at a given time t in the pulse se- 
quence and following (f>2 corresponds to the CID pathway. 
In this case, adiabaticity is achieved if Aqtaf (t) ^ 1 for 
all t. This method has the advantage that an estimate 
for the minimum pulse time required is obtained without 
explicitly integrating the time dependence. The disad- 
vantage is that it provides no information on the fidelity 
(probability of successful transfer) of a given pulse se- 
quence, which still requires numerical integration. Fig. [8] 
shows fmax^CTAp(i) as & fraction of the pulse time for 
each of the viable pulse sequences. Looking at the maxi- 
mum of the adiabaticity, we expect that pulse IV can be 
applied up to 4 times faster than pulse I. It is interest- 
ing to note that the adiabaticity predicts a difference in 
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an electron from well 1 to 3 without populating well 2 
(according to the FS model), we can define the proba- 
bility of deviation from the ideal result as the maximum 
occupancy of the central well and surrounding barriers 
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FIG. 8: ^cTAp(t) plotted for three different pulse sequences 
(pulse III is not included as it requires excessively long pulse 
times). From this plot we can estimate that for the cho- 
sen values of Vmax and Knin, the pulse time must be at least 
tmax > l/f00(i5*)~^ to ensure adiabatic evolution. 



wells 



L/2+W 



L/2-W 



^p*^pdx. (16) 



In Fig. [HI we plot this probability as a function of both 
barrier height Vi = V2 = Vq and barrier width w. 
This plot shows that the population is exponentially sup- 
pressed but non-zero for finite barrier, as foreshadowed 
earlier. 

We can also solve this problem analytically in the limit 
of large barriers. In the limit that the barriers are high or 
wide (or both), the major contribution to P2 comes from 
the exponential decay of the wavefunction as it enters 
the barriers (rather than population in the central well 
itself). We can therefore model this using the solution 
for a particle in a box with one infinite wall and one 
finite wall. In this case, the solution is a transcendental 
equation. 



k2 



tan(fciL) 



(17) 



the behaviour of pulses I and II, despite the similarity 
of the corresponding eigenspectra. We also see that the 
maxima of ^ctap(0 are not at the middle of the pulse 
sequence, in contrast to the FS model. 



CENTRAL WELL POPULATION 



where fci = yj2mE/h?, ^2 = ^2m(Vb — E)IK^ and L is 
the width of the well. The central region probability is 
then approximately the fraction of the particle distribu- 
tion function which is inside the barrier. This approx- 
imate solution is also plotted in Fig. [9] and is valid for 
ui X Vb > lOii?*, i.e. barrier areas greater than about 10 
in normalised units. 



One of the defining features of CTAP protocol in the 
FS formalism is that the central well is unpopulated in 
the adiabatic limit. One key question is whether this fea- 
ture is preserved when including the effects of realistic ID 
potentials. Using the finite state formalism, the popula- 
tion is predicted to be exactly zero for infinitely slow evo- 
lution due to the precise cancellation of the components 
of the wavefunction which correspond to the middle well. 
Therefore the fraction of population in the central well 
due to the finite wells represents an important difference 
between the FS and SW descriptions. When performing 
simulations on the effects of decoherence and fabrication 
defectsS^i^^, it is essential to include these deviations from 
the FS picture as these represent one of the principle er- 
ror modes. Non-adiabatic corrections will also introduce 
some population in higher lying states (as seen in the con- 
ventional Landau-Zener problem) although this effect is 
replicated, at least qualitatively, in the FS formalism. 

Using the exact diagonalisation over a finite lattice, 
we can see that for a finite barrier the population is not 
zero, even in the perfectly adiabatic limit. This non-zero 
population results in the system now being susceptible 
to decoherence and other effects which depend on the 
population of well 2. As the CID pulse should transfer 



VI. TIME EVOLUTION 

To calculate the time evolution of an electron in a 
spatially and temporally varying potential, we use the 
Crank-Nicolson method^i^>^ to approximate the prop- 
agator over a time period At. The time evolution of the 
system is then expressed as 

^'(x,i + At) = C/(Ai)*(x,t), (18) 

where the propagator in Cayley form^i^ is 



2-i{At)H 
2 + i{M)H 



(19) 



This allows the time evolution of the system to be calcu- 
lated for an arbitrary combination of initial state, poten- 
tial and pulse scheme. While this method is much more 
computationally intensive than solving the same prob- 
lem using the FS formalism, it does describe the effects 
of finite barriers, spatially distributed wavefunctions and 
arbitrary potential profiles. 

In Fig. (TU] we plot a series of snapshots showing the 
probability distribution, |4'(x,t)p, as a function of time 
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FIG. 9: Occupation probability of the central region as a func- 
tion of barrier height and barrier width. Note that for suf- 
ficiently high barriers, the occupancy is largely independent 
of barrier width. Even with the barrier at its minimum value 
during a CTAP pulse, the population of the central well and 
the barriers is less than 0.1 for a barrier width of w = L/20. 




FIG. 10: Example trace showing evolution of the probability 
distribution as a function of time during a CID pulse for pulse 
II. Qualitatively similar behaviour is observed for the other 
pulse sequences. 



during the application of pulse sequence II. As the bar- 
riers are lowered, the particle starts to penetrate into the 
barriers (indicated by a spreading of the wavefunction) 
and then ultimately move from well 1 to well 3 with only 
an exponentially suppressed population in the well 2. 

To study the approach to the adiabatic limit, in Fig.fTT] 
we plot the final population distribution |^'(a;, imax)P as 
a function of ^max, using pulse II scaled appropriately for 
each tinax- We identify four regions based on the systems 
qualitative properties, the simplest being regions (a) and 




FIG. 11: Plot showing the final state of the system as a func- 
tion of the total time of pulse II. Note the four key regions 

(a) , (b), (c) and (d) described in the text, which correspond 
to four regions of qualitatively different behaviour. 

(d). In region (a), the pulse sequence is applied so quickly 
that the system doesn't have time to respond, resulting in 
no evolution. In region (d), the system evolution is slow 
enough that the entire transfer is adiabatic, resulting in 
complete transfer from well 1 to 3. In the intermedi- 
ate regions (b) and (c), the transfer is slow enough to 
allow evolution but still non-adiabtic in that the other 
eigenstates are populated during the pulse, resulting in 
imperfect transfer. It is interesting to note that in region 

(b) , all 3 wells are populated at the end of the transfer, 
whereas in region (c) the final state is only composed of 
population in well 1 and 3. This population asymme- 
try arises because the spacing between the eigenstates is 
non-uniform, due to the offset in the central well poten- 
tial. Referring back to Fig. [5l we see that the energy gap 
between the 2nd and 3rd eigenstates is larger than that 
between the 1st and 2nd. We can therefore interpret re- 
gion (b) as the point at which the transfer is fast enough 
to populate all three eigenstates whereas the slower pulse 
in region (c) can only populate the lower two states. As 
the lower two states are only adiabatically connected to 
the 1st and 3rd wells, this results in the observed popu- 
lation distribution. 

To compare the various pulsing schemes, we plot the 
probability of successful transfer as a function of the pulse 
time (tmax) in Fig. 1121 following the similar analysis in 
Rcf. [3. The ordering is as expected from Fig. [5] with a 
pulse time of imax ^ 10"^ allowing adiabatic transfer for 
pulse II whereas pulse IV requires approximately 4 times 
longer. While the fidelity for pulse II & IV are monotoni- 
cally increasing functions, the fidelity function for pulse I 
has a noticeable oscillation. This is consistent with pulse 
I being an approximation to pulse II. Previous work has 
shown that the fidelity function of a CTAP pulse (using 
Gaussian modulation) is strongly dependent on the spac- 
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FIG. 12: Fidelity of the pulse as a function of the total pulse 
time, shown for three different pulses. The FS solution to this 
problem is also shown (solid line) with Gaussian overlapping 
pulses as well as the same problem with a timing mismatch 
(dashed line) such that the time separation between the pulses 
is doubled. The oscillations in the fidelity stem from this 
imperfect overlap in the pulses. 



ing between the pulses^!. In this case the approximation 
used to derive pulse I is only strictly valid near the min- 
imum barrier height (maximum tunnelling rate) whereas 
the overlap between the pulses is strongly controlled by 
the value of the barriers at t = tmax/2. This variation in 
the overlap leads to oscillations in the fidelity function. 

For comparison, the solution of the same problem using 
the FS formalism with appropriate energies is also shown 
in Fig. [T21 For the pulsing scheme given by Eq. 
the fidelity increases monotonically (solid line), although 
the maximum gradient of this line is larger than that 
predicted using the SW formalism. The fidelity obtained 
using the FS formalism is also plotted for a Gaussian 
barrier modulation where the separation between pulses 
has been doubled (dashed line), illustrating the onset of 
oscillations in the fidelity due to mismatch in the pulse 
timing. 

The choice of pulse modulation is therefore important 
with the maximum transfer speed ultimately controlled 
by the adiabatic criteria for the transfer. To optimise 
the transfer rate, it would therefore be useful to design 
a pulse sequence which minimises the adiabatic crite- 
ria while still connecting the initial and final states via 
an adiabatic pathway^^. In practice, this optimisation 
would also need to take into account any restrictions on 
the rate at which the barriers can be varied as well as 
the maximum and minimum values the barrier height 
can take. 



VII. CONCLUSION 

Recent progress in theoretical and experiment control 
of quantum systems has highlighted the link between 
quantum optics and coherent solid-state devices. Many 
of the key developments have benefited from the close 
mathematical and conceptual relationship between these 
fields. This has allowed a fast uptake of ideas, where 
a concept in one field can be quickly mapped to the 
other. We have investigated aspects of this link, specif- 
ically where they apply to the use of adiabatic passage 
techniques in solid-state systems. The characteristics of 
the eigenspectrum and evolution of an electron in a triple- 
well system was studied for different time varying poten- 
tials. These results are directly applicable to electron 
shuttling in confined electron systems such as quantum 
dots or between donor sites. 

While the analogy between the two formalisms stud- 
ied is strong enough to guide future work, care must be 
taken, especially when using realistic potential profiles 
to predict the outcomes of future experiments. For more 
sophisticated and accurate calculations, it will be ulti- 
mately necessary to include more physics than is con- 
tained in the finite state approximation. This is partic- 
ularly evident in the example considered in this paper 
when trying to determine the probability of success for 
various pulse designs and switching times at the limit 
of the adiabatic window. At this point the spatial wave 
model predicts both qualitatively and quantitatively dif- 
ferent behaviour to the finite state model, in terms of 
both transfer fidelity and failure modes. In these situa- 
tions, the solution of the spatially varying wavefunction 
is essential for an accurate description of the system dy- 
namics. 

The goal of our work is to assess the validity of the 
usual (quantum optics motivated) approach to the CTAP 
process where the spatial dependence of the phenomenon 
is ignored. Our model analysis, using the explicit solu- 
tion of the time-dependent Schrodinger equation, demon- 
strates the validity of the usual CTAP approach, point- 
ing out at the same time some of the limitations arising 
from the neglect of the spatial-dependence in the stan- 
dard (and extensively used) CTAP theory. We are partic- 
ularly interested in the application of CTAP in the solid 
state quantum information processing in quantum dot 
or semiconductor donor systems. Since such solid state 
CTAP experiments arc still very far in the future, we have 
used a minimal model so as not to complicate the con- 
ceptual issues by unnecessary numerical details. A real 
experimental situation will necessitate going beyond this 
minimal model, which may turn out to be complicated 
in real solid state systems. 
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